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A STATISTICAL STUDY OF ABILITY OF SIMULATED PILOTS 
TO CONTROL THE APPROACH AND DESCENT OF THE 
LUNAR MODULE TO THE LUNAR SURFACE 

By Jacob H. Lichtenstein 
Langley Research Center 

SUMMARY 

An analytical study was made to assess the ability of a pilot, simulated by permis- 
sible error performance, to control the entire lunar approach and descent with relatively 
simple guidance schemes. The control task consisted of applying retrothrust during the 
hyperbolic approach to establish a lunar orbit. During this lunar coasting orbit, the lunar 
module was separated from the command and service module. A second deorbit thrust 
period was used on the lunar module to establish a coasting descent, and a final thrust 
period was used to put the vehicle into a reasonable landing situation. Throughout the 
maneuver the pilot was permitted to make reasonable errors, selected in a random man- 
ner, in the thrust time and attitude control. The control schemes consisted of main- 
taining a constant vehicle attitude with respect to the line of sight to the lunar horizon 
for the first and second thrust periods, and a constant angle with respect to the line of 
sight to the orbiting command and service module for the third thrust period. 

The results show that a relatively accurate pilot, one whose standard deviation of 
error is 0.1 sec and 0.1° or less, would be necessary to avoid excessively large errors. 
Even such a pilot would be required to switch to another form of terminal guidance to 
bring the flight conditions within acceptable values near touchdown. 

Correction coefficients were determined which permitted corrections to be applied 
to the second and third thrust periods for errors made in the previous thrust periods. 
When these correction coefficients were applied, the final standard error in altitude for 
the worst pilot was reduced from about 25 000 to 3000 ft (7620 to 914 m) and the number 
of impacts with the ground was reduced from 21 to 2 in the total of 51 runs. 

INTRODUCTION 

Simple techniques for guidance, navigation, and control are of interest in space 
flight as a means of reducing hardware complexity or for use in manual backup modes. 
Much work has been done in this area at the NASA Langley Research Center as well as 



in various other organizations ( see refs. 1 to 6, for example) . In most studies of this 
type individual tasks have been considered, either analytically or by simulation. Such 
tasks have included establishment of lunar orbits, lunar landings, rendezvous, and launch 
from the moon. The control tasks in each of these phases have involved orientation of 
the spacecraft and application of a specified thrust level for some time interval. The 
analyses of references 1 to 6 have indicated that simple visual cues, such as views of 
the lunar horizon or of an orbiting spacecraft, could be used for spacecraft orientation 
during the thrusting periods. The analyses have also indicated that the various thrusting 
phases could be accomplished quite accurately and economically. 

The purpose of the present study was to examine analytically how well a pilot, using 
the simplified guidance techniques suggested in the references mentioned previously, can 
perform the sequential tasks required during the lunar landing mission. This mission 
starts with the spacecraft approaching the moon in a hyperbolic orbit and terminates with 
the lunar module close to touchdown. The tasks required of the pilot are to control the 
vehicle thrust in order to (a) establish a near -circular coasting orbit, (b) brake from the 
coasting orbit into a descent transfer orbit, and (c) brake from the transfer orbit to the 
final near -touchdown conditions. This sequence of events is shown in figure 1. Since 
statistical experimental error data for pilots performing the required tasks for the entire 
trajectory are not available, a mathematical model representing the pilots was required. 
This model was defined by the magnitude of the errors in thrust control angle and in the 
response times for the start and stop of thrust. Eighteen pilots were represented by 
various combinations of the magnitudes of the angle and time errors. Statistical data 
were obtained by making 51 complete runs for each pilot, with the actual magnitudes of 
the errors selected by a Monte Carlo technique. 

Since it was recognized that the pilot would make control errors during the descent 
and that the final altitude error may be quite large, some method of compensating for the 
errors may be required. Therefore, a set of error coefficients was developed which 
would indicate the changes to be made during the second and third thrust periods to com- 
pensate for errors made in the previous thrust periods. 

SYMBOLS 


In this investigation, measurements were made in U.S. Customary Units but they 
are also given in the International System of Units (SI). The symbols used in this report 
are defined as follows: 


A 


average value of a parameter, 



n 


C correct value of a parameter obtained from nominal rim 
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h 

l 

M 


P 

R 

R o 

S 

t 

*0 

*1 

*2 

*3 

*4 

*5 

*6 

*7 


altitude above the lunar surface, ft (m) 
distance between the LM and the CSM, ft (m) 
mass of the moon, slugs (kg) 
mass of the vehicle, slugs (kg) 
initial mass of the vehicle, slugs (kg) 

number of measurements of the parameter under consideration 
value of the parameter under consideration 

radius from the center of the moon to the vehicle, n. mi. or ft (m) 
radius of the moon, 5 702 000 ft (1 737 969.6 m) 

standard error of a parameter, 
time, sec 

time at start of computations (at h = 200 n. mi. or 370 400 m) 
time at start of first thrust period, sec 
time at end of first thrust period, sec 
time at start of second thrust period, sec 
time at end of second thrust period, sec 
time at start of third thrust period, sec 
time at which the vehicle reaches minimum altitude, sec 
time at which the vehicle reaches the moon's surface, sec 
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Ati 


Atn 


At m 

v 

Xi,Yi 


x i>Yi 


X,Y,Z 


a 


y 


^1 


^2 

e 

\ 

k 


a 


T 


error in the length of the first thrust period, sec 
error in the length of the second thrust period, sec 
time error at the start of the third thrust period, sec 
velocity of the vehicle, ft/sec (m/sec) 
inertial axes (see fig. A- 1(b)) 

distances along and axes, respectively, ft (m) 

body axes of the vehicles, oriented as shown in figure 2, with the subscript 
denoting the vehicle 

angle between the local vertical and the line of sight from the LM to the 
CSM, deg 

flight -path angle, measured up from the local horizontal, deg 

angle between the line of sight from the CSM to the trailing lunar horizon 
and the X^g^j axis during first thrust, deg 

angle between the line of sight from the LM to the trailing lunar horizon 
and the Xl M axis during second thrust, deg 

i y* 

angular displacement of the vehicle in the X;Y; plane, tan -1 -1 rad 

i. *, 

longitude on the lunar surface, deg 

3 

universal gravitational constant, 6.668 x 10”^ km 

g sec' 5 

angle between the X body axis and the local horizontal, deg 

\jjip - a ) 2 

standard deviation of a parameter, U — — — j — 

angle between the local vertical and the line of sight to the lunar horizon, deg 
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V 


angle between the line of sight from the LM to the CSM and the Xj^ axis, 
deg 


Subscripts: 


A 

altitude 

a 

control angle 

CSM 

command and service module 

LM 

lunar module 

Roe 

range 

t 

time 

V 

velocity 


Dots above a symbol denote differentiation with respect to time. 

METHOD 

The analysis and results of this paper are based on the computed ability of pilots 
to fly along a predetermined nominal trajectory by adhering to a time-defined nominal 
flight plan. It is assumed that the spacecraft is approaching the moon along a hyperbolic 
trajectory similar to that for the Apollo mission, and the nominal trajectory starts at an 
altitude of 200 n. mi. (370 400 m) above the surface. The terminal conditions are an 
altitude of 1 n. mi. (1852 m), zero vertical velocity, and a forward velocity of about 
600 ft/sec (182.9 m/sec). In order to obtain data for the statistical analysis, 51 descent 
trajectories in which the control errors were randomly selected were computed for each 
of the 18 pilots considered. 

The trajectory computations were made by using a general three -degree -of -freedom 
system of equations with the moon as the central body. For this investigation, however, 
it was decided that there would be no loss in generality if the problem were restricted to 
the planar case with the vehicle free to move longitudinally and vertically and to rotate 
for thrust alinement. The computations permit controlled thrust of the combined vehicle 
(command and service module with attached lunar module) during its approach, separation 
of the two vehicles, and thrust of the lunar module while keeping track of the orbiting 
command and service module. Because this was a general investigation, no attempt was 
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made to duplicate a specific Apollo mission, in which the lunar horizon may not be visible. 
For this investigation, it was assumed that the lunar horizon would be visible and that 
there would be no uncertainty in horizon determination due to terrain irregularities. A 
more detailed discussion of the equations used during the pilot control part of the experi- 
ment is presented in appendix A. A separate computational program was used to reduce 
the resulting statistical data and to make the plots of these data. 

Nominal Flight Path 

The nominal flight path was determined by an iteration procedure which resulted in 
the desired trajectory characteristics at key phases of the mission. This trajectory was 
designed to allow manual control of the approach, and no attempt was made to duplicate 
the present Apollo approach. The nominal trajectory starts at 200 n. mi. (370 400 m) 
along the hyperbolic approach, and at that point it has the following parameters: 

V = 7945.04 ft/sec (2421.65 m/sec) 

y = -20.722 936 5° 

h = 1 215 406.2 ft (370 455.8 m) 

Three thrust periods are used to accomplish the mission. The first slows the combined 
vehicle from the approach velocity and establishes a nearly circular orbit at 80 n. mi. 

(148 160 m). The second slows the lunar module from orbital velocity to establish a 
descent transfer orbit with periapsis at 50 000 ft (15 240 m). The third is the final slow- 
down to the terminal conditions. These thrusting tasks are depicted in figure 1. The 
body axes of the vehicles are oriented as shown in figure 2. 


Start of problem 
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(a) Establishment of circular orbit. 

Figure 1.- The three thrusting tasks required of the pilot. 




Command and service 
module 



(b) Establishment of transfer orbit. 


(c) Landing approach. 


Figure 1.- Concluded. 



Figure 2.- Sketch of the command and service module and the lunar module coupled 
together, with positive directions of their body axes indicated by arrows. 



In order to follow this trajectory the pilot's mission time schedule is as follows: 


tg = 0.000 sec 

Problem is initiated at 200 n. mi. (370 400 m). Remain in the coasting orbit, 
tj = 205.000 sec 

Aline the vehicle so that the X^g M axis is at an angle of 50.9° with the receding 
lunar horizon (see fig. 1(a)). Initiate thrust and maintain until t 2 - 

t2 = 536.100 sec 

End thrust. Vehicle should now be in a nearly circular coasting orbit at 80 n. mi. 
(148 160 m). 

t 2 to t 3 

Coast in nearly circular orbit, during which separation of LM and CSM occurs. 
t 3 = 6580.000 sec 

Aline the LM so that the X LM axis is at an angle of 67.15° with the receding lunar 
horizon (see fig. 1(b)), and initiate thrust in the LM. This time was chosen because 
it is at the periapsis of the nearly circular orbit. 

t 4 = 6589.445 sec 

End thrust. The LM should be in a descending transfer orbit. 
t 4 to t 5 

Coast in transfer orbit 
t 5 = 10 069.000 sec 

Aline LM so that the X LM axis is at an angle of -14.75° with the line of sight to 
the orbiting CSM and start thrust (see fig. 1(c)). Maintain this angle throughout the 
thrust period. This time was chosen because it is the periapsis of the transfer orbit 
at h = 50 761 ft (15 472 m). 


8 


tg = 10 462 sec 

End thrust. The LM should be at the terminal conditions: 

h = 6465 ft (1970.5 m) 
h *0 

V = 549.3 ft/sec (167.4 m/sec) 

These end conditions are similar to those for a medium -speed airplane in level flight, 
and it is assumed that a pilot with proper training can make a safe landing on the lunar 
surface. 


Pilot Performance Definition 

Since this was an analytical study, a definition for the mathematical model of the 
pilot was required. For any run (single trajectory computation) the magnitude and direc- 
tion of the possible errors would occur in a random manner. The spectrum of the errors 
for any task is defined by the error distribution and standard deviation. In this investiga- 
tion it was assumed that the pilot's errors would be normally distributed about zero. The 
standard deviations (cr a and o^) were based on some unpublished tests using an engineer 
as a pilot. These tests were multitask jobs in which the pilot's attention was applied pri- 
marily to maintaining three -axis attitude control, and initiation and cutoff of thrust were 
in response to an audible signal. The results showed that starting and stopping the thrust 
could readily be controlled to 0.1 sec and the angle could be controlled to within 0.2°. A 
total of 18 pilot definitions were used in this study, and they are specified in table I in 
terms of a a and oj.. 


TABLE I.- PILOT PERFORMANCE PARAMETERS 


Pilot 

<x a , deg 

Oj., sec 

1 

0 

0.05 

2 

0 

.10 

3 

0 

.20 

4 

.05 

0 

5 

.05 

.05 

6 

.05 

.10 

7 

.10 

0 

8 

.10 

.05 

9 

.10 

.10 

10 

.10 

.15 

11 

.10 

.20 

12 

.15 

.10 

13 

.15 

.15 

14 

.15 

.20 

15 

.20 

0 

16 

.20 

.10 

17 

.20 

.15 

18 

.20 

.20 
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A random -number computing subroutine, called the Monte Carlo technique, was used 
to obtain the errors at each of the actions required of the pilot. The actual distribution 
of the errors as they occurred during a typical series of runs (cr a = 0.20°, cr^ = 0.20 sec) 
is shown in figure 3. In figure 3(a) the number of occurrences of the errors in each mag- 
nitude bracket is shown for each time (tj, t£, etc.) and for the compilation of all five 
times. Similar data for the angle errors are shown in figure 3(b). Because of the rela- 
tively low number of samples (51 runs) at any time point, the error distribution is only 
approximately normal; however, if all the timing errors for a series (errors at tj, t2, 
t3, 1 4, and t5 for all 51 rims) are considered, it can be seen that the error distribution 
is very close to a normal distribution. The nearly straight line that results when the 
cumulative percentage total is plotted on probability paper confirms this conclusion. 


Vehicle 

The vehicles assumed in this investigation were similar to vehicles considered for 
the Apollo mission. A sketch of the command and service module and the lunar module 
is shown in figure 2. The mass and thrust magnitude data for the combined vehicle 
(command and service module with the attached lunar module) were obtained from the 
Apollo literature available at the time this investigation was started. The values used 
in this investigation for the combined and lunar modules are: 


Combined vehicle: 

Mass 

Thrust 

Mass flow during thrust 


. 2741.67 slugs (40 011.7 kg) 

. . . . 21 897 lb (97 402.7 N) 
2.259 slugs/sec (33.0 kg/sec) 


Lunar module: 

Mass 974.13 slugs (14 212.6 kg) 

Thrust 10 000 lb (44.482.2 N) 

Mass flow during thrust 1.0318 slugs/sec (15. 1 kg/sec) 


The direction of thrust is along the X^g M axis in the positive direction for the com- 
bined vehicle and along the axis in the negative direction for the lunar module. 

Since this investigation was concerned with piloting performance, errors such as thrust 
misalinement, flow irregularities, and thrust tailoff were not considered. 
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Number of occurrences 




-.50-45-40 r35 -30:25 - 20 - 1 5 -I0-.05 0 .05 .10 .15.20.25 .30.35 4 0 45.50 .55 .60 


Magnitude of error, deg 

(b) Distribution of the angle errors. 
Figure 3.- Concluded. 





RESULTS AND DISCUSSION 


Standard Flight Procedure 

Each computed trajectory was examined at several points to determine departures 
in altitude, velocity, and range from the nominal trajectory. The points examined were 
the end of the first and second thrust periods ^2 and t 4 , respectively) and the start and 
end of the third thrust period (t§ and tg, respectively). The third thrust period was 
assumed to terminate either when the rate of descent became zero or on ground contact, 
whichever occurred first. Examples of the type of data obtained are given in appendix B. 
These data are in the form of bar graphs showing the distribution of errors in velocity 
and altitude for the various pilot definitions. The results were compiled in figures 4 to 8 
in a form more suitable for analysis. These figures show the effects of pilot definition 
(a a and cr t ) on the magnitude of the stardard errors in altitude, velocity, and range. 






(a) Pilot definition angle error. 

Figure 4.- Variation of the standard error in altitude with pilot definition. 


13 


-]280x|tf 





iff' 

Hi 60 g 


STD OEV OF PILOT DEFINITION TIME. CT t . SEC 



-.1 100x10' 
-1000 
- 900 


A 800 


e 


700 - 

600 g 

500 J 
< 

fe 

400 

§ 

a 

- 300 u 
a 

■ 200 


id 00 


(t>) Pilot definition time error. 

Figure 4.- Concluded. 

Altitude errors . - The standard error in altitude as a function of pilot definition is 
shown in figure 4 for several key times along the trajectory. Since all the trajectories 
start with the same initial conditions and the standard deviations of the pilot errors (a a 
and erf.) are reasonably small, there is only a modest buildup of the error up to the second 
thrust period. At this time the largest error was about 1 n. mi. (1852 m) in the 80 n. mi. 
(148 160 m) coasting orbit. However, this error combined with the error in velocity at 
the end of the second retrothrust can cause the largest buildup of errors to occur during 
the descending transfer orbit. These error buildups are large enough to jeopardize the 
landing unless corrective action is taken. The altitude errors presented for time tg 
are for the trajectories that did not hit the lunar surface and are errors in the minimum 
altitudes achieved during this thrust period (rate of descent is zero at tg). These errors 
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are essentially positive, as large negative errors would cause the vehicle to hit the sur- 
face. The seriousness of large negative errors is indicated in figure 5, where each point 
represents the number of surface impacts out of the 51 runs of a pilot. The data show 
that 270 of the total of 918 runs, or about 30 percent, hit the ground. Even for relatively 
small pilot errors (standard deviations of 0.10 or less) there is still about a 20-percent 
chance of hitting the surface if no correction is made to the flight procedure. The data 
for the standard error of the vertical velocities at impact, presented in figure 6, show 
that these velocities approach 300 ft /sec (91.4 m/sec). 



Std. dev.of pilot definition angle, a Q , deg 


(a) Pilot definition angle error. 

Figure 5 • - Variation of the number of times the lunar module hit the ground with variation 
in the pilot definition. A total of 51 runs was made for each pilot. 
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Standard dev. of pilot definition angle, deg 

O cr a = 0 

□ <?n - .05 

O oq = .10 

A Oq = .15 

CTq = .20 

oq = . 20 with corrections 




Std, dev. of pilot definition time , 


(t>) Pilot definition time error 


Figure 5.- Concluded. 



STANDARD ERROR OF VERTICAL VELOCITY AT IMPACT, FT/SEC 



Figure 6.- Vertical velocity at the time of impact with the ground as a 

function of pilot definition. 
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STANDARD ERROR OF VERTICAL VELOCITY AT IMPACT, m/SEC 



Velocity errors. - The buildup of the standard error of the velocity as the trajec- 
tories progress is shown in figure 7. The velocity errors are reasonably small even up 
to the start of landing retrothrust at tg, where the error is only 14 ft/sec (4.3 m/sec) 
at a velocity of 5582 ft/sec (1701.4 m/sec). It is during this last thrust period (between 
ts and tg) that the largest errors are developed. Note that the data presented at tg 
are for the runs that did not impact, and are thus at minimum altitudes where h = 0. 

The errors, therefore, are longitudinal velocity errors. 




STD DEV OF PILOT DEFINITION TIME. SEC 
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(a) Pilot definition angle error. 

Figure 7*- Variation of the standard error in velocity with pilot definition. 
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Range errors . - The standard errors in range are shown in figure 8. There is not 
much buildup in range error as time progresses. Most of the error is acquired during 
the first thrusting period, and generally is less than 5 n. mi. (9260 m). 




(a) Pilot definition angle error. 




(b) Pilot definition time error. 

Figure 8.- Variation of standard error in range with pilot definition. 
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Correction of Terminal Altitude 


Probably the most important conditions to be met in this study are the terminal 
altitude of about 6500 ft (2000 m) with zero rate of descent. The results of this study 
indicate the possibility of appreciable departure from these conditions by the assumed 
pilots following the specified procedure. It is therefore important to obtain, early in the 
task, some estimate of what the terminal conditions will be, and to correct the trajectory 
to obtain the desired terminal conditions. This information can be obtained during the 
relatively long coast times after the first and second thrust periods. 

An analytical study was made to determine the effects on the terminal conditions of 
thrust aiming errors and timing errors for each of the thrust periods individually while 
the other thrust periods adhered to the nominal plan. The results, presented in figure 9, 
show that the terminal altitude error is nearly a linear function of the individual errors 
in angle or time. As a result of this near linearity, the results can be expressed as sen- 
sitivity coefficients as follows: 


1. Final altitude error per degree of first -thrust angle error (fig. 9(a)): 

= -92 500 ft/deg (-28 194 m/deg) 

d(A Vl ) 

2. Final altitude error per second of first -thrust time error (fig. 9(b)): 


■^4= -38 100 ft/sec (-11 612.9 m/sec) 

8 ( At l) 

3. Final altitude error per degree of second-thrust angle error (fig. 9(c)): 
= -3020 ft/deg (-920.5 m/deg) 


a(Ah) 


4. Final altitude error per second of second-thrust time error (fig. 9(d)): 


a(Ah) _ 

8 ( At n) 


-65 900 ft /sec (-20 086.3 m/sec) 


5. Final altitude error per degree of third-thrust angle error (fig. 9(e)): 


8(Ah) 

d(Av) 


= +14 200 ft /deg (4328.2 m/deg) 


6. Final altitude error per second of time error in starting third thrust: 
= 251.4 ft/sec (76.6 m/sec) 

9( At m ) 
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Figure 9 * - Continued. 




These control coefficients offer a means of rectifying errors made during a thrust 
period by applying corrective measures in the subsequent thrust periods. If an angle or 
time error was made during a thrust maneuver, the coefficients could be used to predict 
the magnitude of the altitude error that would result at the end of the trajectory if the 
nominal plan is followed. If the magnitude of the error is unacceptable, corrective mea- 
sures are necessary. By examination of the coefficients for the following thrust periods, 
the corrections to the angle or time which would be required to compensate for the initial 
error can be determined. Exact correction of the altitude error probably will not be 
obtained because the curves in figure 9 are not exactly linear and because the error 
coefficients were determined as sensitivity parameters about the nominal trajectory, 
whereas they will be applied to off-nominal conditions. 



Figure 10.- Residual altitude errors after correction of first-thrust time errors 
hy corrections during second and third thrust periods. 


The results of some computations using these coefficients to correct for definite 
time and angle errors imposed during the first thrust period are presented in figures 10 
and 11, and summarized in table n. These show the improvement that can be obtained 
for an individual error. 

In order to evaluate the benefit to be obtained by using these corrections in trajec- 
tories where the errors are random in magnitude and can be made throughout the trajec- 
tory, these correction coefficients were inserted into the basic computation. Pilot 18 
(°a = °*20°; cr t = 0.20 sec) was chosen as the example, and the correction scheme used 
was as follows: 
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1. For both angle and time errors made during the first thrust, a new nominal 
second -thrust time was computed by adjusting the cutoff time. 

2. For both angle and time errors made during the second thrust, a new nominal 
third -thrust angle was computed. 

3. Errors incurred during the third thrust remained uncorrected since there was 
no subsequent thrust period in which corrections could be made. 



Figure 11.- Residual altitude errors after correction of first -thrust angle 
error by corrections during second and third thrust periods. 
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TABLE H.- EFFECTIVENESS OF COMPENSATING FOR TERMINAL ALTITUDE ERRORS 


Procedural 

error 

Error in h due to 
procedural error 

Compensating 

adjustment 

Estimate of error 
in corrected 
terminal altitude 

Atp 

A Jv 

feet 

meters 

At n» 

t”2’ 

Ap, 

feet 

meters 

sec 

deg 



sec 

deg 

deg 



+0.4 

0 

-15 240 

-4 645.1 

0 

-5.06 

0 

2 200 

670.6 

+.2 

0 

-7 620 

-2 322.6 

0 

-2.53 

0 

550 

167.6 

-.2 

0 

7 620 

2 322.6 

0 

+2.53 

0 

100 

30.5 

-.4 

0 

15 240 

4 645.1 

0 

+5.06 

0 

300 

91.4 

+0.4 

0 

-15 240 

-4 645.1 

-0.232 

0 

0 

-200 

-61.0 

-.4 

0 

15 240 

4 645.1 

+.232 

0 

0 

-100 

-30.5 

+•4 

0 

-15 240 

-4 645.1 

0 

0 

-1.08 

-400 

-121.9 

-.4 

0 

15 240 

4 645.1 

0 

0 

+1.08 

-300 

-91.4 

0 

+0.4 

-37 000 

-11 277.6 

0 

-12.28 

0 

11 200 

3413.8 

0 

+.2 

-18 500 

-5 638.8 

o 

-6.14 

0 

2 900 

883.9 

0 

-.2 

18 500 

5 638.8 

0 

+6.14 

0 

3 400 

1036.3 

0 

-.4 

37 000 

11 277.6 

0 

+12.28 

0 

14 300 

4358.6 

0 

+0.4 

-37 000 

-11 277.6 

-0.563 

0 

0 

300 

91.4 

0 

-.4 

37 000 

11 277.6 

+.563 

0 

0 

500 

152.4 

0 

+.4 

-37 000 

-11 277.6 

0 

0 

+2.61 

-950 

-289.6 

0 

-.4 

37 000 

11 277.6 

0 

0 

-2.61 

-400 

-91.4 
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The results of these computations are shown in figures 4 to 8 as the solid symbols 
at cr a = 0.20° and = 0.20 sec. The improvement in the final altitude is shown in 
figures 4 and 5; the standard error at time tg decreased from about 25 000 ft (7620 m) 
to 3000 ft (914.4 m), and the number of times the vehicle hit the ground decreased from 
21 to 2 in the total of 51 runs. The vertical velocity at impact also decreased signifi- 
cantly, from about 240 ft/sec (73.2 m/sec) to less than 160 ft/sec (48.8 m/sec). These 
residual errors result mainly from the control errors incurred during the third thrust 
period. The other end conditions - total velocity and range errors - were not signifi- 
cantly affected. 

The results presented in figures 10 and 11 and table II show that terminal altitude 
errors can be substantially decreased by using the simple error -coefficient method of 
compensation. In addition, the results show that corrections to the second thrust time or 
the third thrust angle are more effective than corrections to the second thrust angle. 

CONCLUDING REMARKS 

An analytical study was made to assess the ability of a pilot, simulated by permis- 
sible error performance, to control the entire lunar approach and descent with relatively 
simple guidance schemes. The control task consisted of applying retrothrust during the 
hyperbolic approach to establish a lunar orbit. During this lunar coasting orbit, the lunar 
module was separated from the command and service module. A second deorbit thrust 
period was used on the lunar module to establish a coasting descent, and a final thrust 
period was used to put the vehicle into a reasonable landing situation. Throughout the 
maneuver the pilot was permitted to make reasonable errors, selected in a random man- 
ner, in the thrust time and attitude control. The control schemes consisted of maintaining 
a constant vehicle attitude with respect to the line of sight to the lunar horizon for the first 
and second thrust periods, and a constant angle with respect to the line of sight to the 
orbiting command and service module for the third thrust period. 

The results show that a reasonably accurate pilot, one whose standard deviation in 
the time and angle control is within 0. 1 second or degree, was required so that the final 
altitude and velocity errors would not become too large. Even with a pilot of these spec- 
ifications, the errors were so large that the vehicle hit the ground about 20 percent of the 
time, with a vertical velocity as high as 260 ft/sec (79.2 m/sec), if no correction was 
made to the flight path. 

A set of coefficients that permit correction of the final altitude error was developed. 
In order to use these coefficients, it is necessary to obtain a measure of the magnitude of 
the errors made during a thrust period so that a proper correction can be applied to the 
following thrust period. As an example of the corrections that could be made by means of 
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these coefficients, suppose that an error of 0.4° was made in the first thrust angle. This 
would incur a penalty of ±37 000 ft (11 277.6 m) in the final altitude; however, corrective 
measures applied during the second and third thrust periods could reduce this error to 
less than 1000 ft (304.8 m). 

When these correction coefficients were used in the basic trajectory computation 
program, the final altitude errors were significantly decreased. For example, for the 
worst pilot (standard deviation of 0.20° and 0.20 sec) the final standard error of altitude 
was reduced from 25 000 to 3000 ft (7620 to 914 m) and the number of surface impacts 
was reduced from 21 to 2 in 51 rims. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 14, 1971. 
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APPENDIX A 


MATHEMATICAL EQUATIONS AND METHODS USED 
IN THE INVESTIGATION 

For this investigation all the orbits were restricted to one plane (see fig. A-l) and 
the vehicle was restricted to essentially three degrees of freedom, motion longitudinally 
and vertically and rotation for thrust alinement. This last condition was not really a 
degree of freedom because the thrust vector was positioned without resort to rotational 
dynamics. 

The equations of motion used in this investigation were 
mxi = Fxi - mG cos X 


myi = Fyi - mG sin X 


where Fxi and Fy^ are the thrust components along the and axes, respec- 
tively. The direction of thrust was different for the two vehicles, being along the Xqsm 
axis for the CSM and along the -Z LM axis for the LM. The expressions for the thrust 
are then: 

For the CSM, 

F Xi = F CSM sin ^ " x ) 

F Yi = F CSM cos ^ " 

For the LM, 

F Xi = F LM cos ^ “ X ) 

F Yi = F LM sin ^ " X ) 

The terms F^g M and Fj^ are considered to be constant. 

The term G is the gravitational vector and is given by 



APPENDIX A - Continued 


The instantaneous mass is given by 
m = m Q - Jmdt 

where m is integrated during the thrust periods. The distance between the vehicles 
was obtained from 



The radial distance from the center of the moon to either of the vehicles is 
R = l/xi 2 +Yi 2 


and the velocity of the vehicle is 

V = j/xj 2 +^2 

Two schemes for controlling the angle £ were used: 

1. Control scheme 1 takes its cue from the lunar horizon and maintains the control 
angle at a constant value with respect to the trailing horizon. The angle £ varies as 
the angle to the horizon varies. The angle r between the local vertical and the horizon 
is shown in figure A- 1(a) and was obtained from 

r = tan"* — — 

h(2R 0+ h) 1/2 

The angle from the X body axis, which is the one that the pilot uses as a refer- 
ence, to the horizon is given by 

77 = -(90° - £ - r) 


and 


£ = ij + 90° - r 

The angle rj is the angle that the pilot maintains at the specified value. 
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APPENDIX A - Concluded 


2. Control scheme 2 takes its cue from the orbiting command module. The geom- 
etry for this control scheme is shown in figure A- 1(b). The angle between the line of 
sight from the LM to the CSM and the local vertical is a, which is obtained from 

cos a = R LM ~ R CSM cos ** 

l 

The control angle v, the angle between the axis of the LM and the line of 

sight to the CSM, is obtained from 

v = a - 90° - £ 

and the thrust control angle is thus 
| = Q!- 90° - v 

The angle v is the angle that the pilot maintains at the specified value. 
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APPENDIX B 

DATA FOR THE ERROR DISTRIBUTIONS 

The basic data of this study were obtained in terms of errors in the trajectory 
parameters at key times along the trajectory. The data presented in figures B-l to B-5 
show the distributions of the errors in altitude, velocity, and range. For pilot 18 (the 
worst one) the distributions are shown for several times along the trajectory (figs. B-l 
and B-2). For the other pilots, however, the data are presented only for the end of the 
run, either the time of minimum altitude or the time that the vehicle hit the ground 
(figs. B-3 and B-5). 

The altitude and velocity error distributions as shown in figures B-l and B-2 for 
pilot 18 are at the end of the first and second thrust periods (t 2 and 1 4 , respectively), 
at the start of the third thrust ^ 5 ) and at the end of the run (tg or t 7 ). Near the start 
of the trajectory (t 2 ) the error distribution is reasonably close to normal, but as the tra- 
jectory proceeds the distribution becomes more nearly uniform except for the number of 
times the LM hits the ground. The maximum magnitude of the altitude errors increases 
very rapidly, progressing from about ±4000 ft (1219.2 m) at t 2 to ±18 000 ft (5486.5 m) 
at t 4 , ±45 000 ft (13 716 m) at tg, and +56 000 ft (17 068.8 m) at tg. The nominal final 
altitude was 6465 ft (1970.5 m); therefore, any error in the negative direction larger than 
6465 ft (1970.5 m) would cause the vehicle to hit the ground. This puts a lower limit on 
the possible altitude. The column on the left end of the scale is so high because it con- 
tains not only the errors that would normally fall in this band, but also all larger negative 
errors. The pattern of the velocity error distributions is similar to that of the altitude 
error distributions. The velocity errors are as large as 6000 ft/sec (1829 m/sec). 

The distributions of the errors in altitude, velocity, and range at the end of the tra- 
jectories (times tg and t 7 ) are shown in figures B-3 to B-5 for the various pilot defini- 
tions. The data are presented to show the effect of changing 0 ^ for a constant magnitude 
of cr a . As expected, the magnitude of the errors increases with increases in both types 
of permissible pilot errors. Of course, as the range of the magnitude of the errors 
increases and spreads out along the abscissa, the number of occurrences for each error 
division decreases; therefore, the error distribution gradually changes from one which is 
nearly normal to one that approaches a uniform distribution across the spectrum. 

The error -distribution data for pilot 18 (cr a = 0.20°, a t = 0.20 sec) with the fir st- 
and second-thrust errors automatically corrected are shown in figure B- 6 . A comparison 
of this figure with figure B-l(d) shows the considerable decrease in the altitude -error 
range that can be achieved with such corrections. 
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APPENDIX B — Continued 



MAGNITUDE OF ERROR. FEET 


1372 -1219 -1067 -914 -762 -610 -457 -305 -152 0 152 305 457 610 762 914 1067 1219 1372 024 1676XI0 1 

MAGNITUDE OF ERROR, METERS 


[c) At the start of third thrust (t^j . 


MAGNITUDE OF ERROR. FEET 

I I I I . . -J I I !_ -I 1 . I I I [ I — l l i - 1 1 I i 

213 -107 6 107 213 320 427 533 640 747 853 960 1067 1173 1280 1387 1494 1600 1707 1814 l920XfO* 

MAGNITUDE OF ERROR, METERS 

(d) At the end of the run (tg or t^) . The large column at the left includes the 

trajectories that impact the ground. 


Figure B-l.- Concluded. 



NUMBER OF OCCURRENCES 


APPENDIX B - Continued 




MAGNITUDE OF ERROR, m/SEC 


(b) At the end of second thrust 


(H)* 


Figure B-2.- Velocity- error distributions for pilot 18 at various times 

along the trajectory. 
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-21 -12 


I 1 1 L 


I I I I I I 1 


1111 1-1 _l I I 1 I I I 


15 21 34 43 52 61 70 73 SB 30 107 116 125 134 143 1 

MAGNITUDE OF ERROR, m 


I I I I I I 1 1 

-21 -13 -16 -13 -10 -0 -5 -2 .61 3 6 3 12 14 17 20 23 25 28 31 

MAGNITUDE OF ERROR, m 



I J I 1 I 

-2t -15 -9 -3 3 


J l l I .1 I I J 

3 15 21 27 34 40 46 52 

MAGNITUDE OF ERROR, m 


64 70 76 82 


J I I I I 


(a) Angle error definition cr a = 0°. 

Figure B-3.- Altitude -error distribution at the end of the run (tg or ty) for the 

various pilots. The large column at the left side includes the trajectories that 
impact the surface. 


X 10* 


x 10* 
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NUMBER OF OCCURRENCES ^ NUMBER OF OCCURRENCES 


APPENDIX B — Continued 



-21 -19 -16 -14 -12 -3 -7 -4 -2 .61 3 S 

MAGNITUDE OF ERROR. 


10 13 15 10 20 23 25 27 x 10* 


-21 -15 -9 -3 3 9 15 21 27 34 40 46 S2 50 

MAGNITUDE OF ERROR, m 



MRGNITUOE OF ERROR. FT 


I I 1 L_L_LJ._IJ _L_J I I- I 1 J_l f 

-21 -15 - 9 - 3 3 9 15 21 27 34 40 46 52 50 64 70 76 82 

MAGNITUDE OF ERROR, m 


Ob) Angle error definition a a = 0.05°. 
Figure B-J.- Continued. 
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NUMBER OF OCCURRENCES ^ , _ NU _ nBER ° F OCCURRENCES ^ NUMBER OF OCCURRENCES 


APPENDIX B - Continued 



MAGNITUDE OF ERROR. FT MAGNITUDE OF ERROR. FT 


L I. I ! I I I I I I I I. J I I I I L_ 

-21 -IS -3 -3 3 9 IS 21 27 34 40 46 52 SB 64 70 76 82 

MAGNITUDE OF ERROR, m 
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NUMBER OF OCCURRENCES NUMBER OF OCCURRENCES 


APPENDIX B - Continued 




MflGNI TUDE OF ERROR, FT MAGNITUDE OF ERROR. FT 

l I I I I l I I l 1 I I I I I 1 I 1 111 1 L I I 1 I 1 I I I 1 J I ! 1 

-2\ -12 -3 6 IS 24 34 43 52 61 70 73 88 38 107 116 12S 134 143 152 162 * 10* -21 ^11 0 11 21 32 43 S3 64 75 05 96 107 117 128 139 143 160 171 101 192 x 10* 

MAGNITUDE OF ERROR, m MAGNITUDE OF ERROR . m 


(e) Angle error definition cr a = 0.20°. 
Figure B-3.- Concluded. 
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NUMBER OF OCCURRENCES NUMBER OF OCCURRENCES 


APPENDIX B - Continued 



20 32 44 SB 63 81 33 105 117 130 142 1S4 166 178 191 203 215 227 

MAGNITUDE OF ERROR . "/SEC 


-2 4 10 16 22 2 6 34 41 47 53 59 6S 71 77 83 89 3S 101 108 114 12Q x 10 l 

MAGNITUDE OF ERROR. m/SEC 



MAGNITUDE OF ERROR. "VSEC 

(a) Angle error definition cr a =0°. 

Figure Velocity-error distributions at the end of the run (tg or tjj for 

the various pilots. 
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I 


APPENDIX B - Continued 




PILOT 5, o} = .05 SEC 


MRGNI TUOE OF ERROR. m/SEC 


72 70 04 30 36 102 100 114 : 

MAGNITUDE OF ERROR. FT/SEC 

— J L. 1 1 .1 t l 1 I l . J J l .1 J I I 

S 7 3 II 13 IS 16 10 20 22 24 26 27 23 31 33 3S ; 

MAGNITUDE OF ERROR. "/SEC 



J I I 1 1 ! L 


■ J I I J 


L L I J__J I I 

-2 .30 2 5 7 9 11 13 15 17 20 22 24 26 2B 30 32 34 37 39 41 

MAGNITUDE OF ERROR. "/SEC 


10' 

10‘ 


(t)) Angle error definition a a = 0.05°. 
Figure B-4.- Continued. 
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Continued 




i i i 

-.61 2 4 


j j .i.j. i j i j i j 

10 12 14 16 19 21 23 25 27 29 31 34 


MAGNITUDE OF ERROR. ™/SEC 


J 

36 


i ] 


J 

42 * 10* 


I I I I 

-2 -91 4 6 


_L 

9 


_L 

12 


-I I J l_.J .11 1 I 

IS 17 20 23 26 20 31 34 37 

MAGNITUDE OF ERROR. m/SEC 


42 45 '48 50 S3 x 10» 




I I I I I I t I I I I I I I I I I I I I I 

-S -2 2 S 0 11 14 17 20 23 26 29 32 35 30 41 44 47 SO 53 56 X 10' 

MAGNITUDE OF ERROR. m/sEC 


1 I L L.J 1 I t I I 1 1 I I 1 I I 1 L I I 

-3 3 9 IS 2t 27 34 40 46 52 58 64 70 76 82 08 94 101 107 113 119 X 10 l 

MAGNITUDE OF ERROR. m/SEC 



-3 3 9 15 21 27 34 40 4 6 52 SB 64 70 76 02 88 94 101 107 113 119 X 10' 

MAGNITUDE OF ERROR. "VSEC 

(c) Angle error definition a a = 0.10°. 
Figure B-4.- Continued. 
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APPENDIX B 


Continued 



PILOT 12, = .IOSEC 


\ .1 

B5 205 225 245 265 

MAGNITUDE OF ERROR. FT/SEC 


I - 1 - I I I — I 


PILOT 13, 0} = 


Tba: 


irOzzl _L._Jr-L 


j i j j i i i i j 


jj I i 


L.J I I I I J . . . _ 

-5 2 8 14 20 26 32 3B 44 50 SB 62 69 75 B1 07 93 99 105 111 117 * 10‘ 

MAGNITUDE OF ERROR. ™/SEC 


^10 10 30 “SO 70 90" HO 130 150 170 190 210 230 250 270 t 
MAGNITUDE OF ERROR. FT /SEC 

UL1.J.J J._L.LJ_IJ J 1 J _1 

-3 3 9 15 21 27 34 40 4B 52 58 64 70 76 82 

MAGNITUDE OF ERROR. "'SEC 



MAGNITUDE OF ERROR. m/SEC 


10' 


io> 


(d) Angle error definition a a = 0.15°. 
Figure B-4.- Continued. 


.15 SEC 


J I I f I I 

»30 310 330 350 370 390 t 10 l 

J J _1_1 I I 

B9 94 101 107 113 119 * 10 l 
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NUMBER OF OCCURRENCES 


APPENDIX B — Continued 


PILOT 15, oj = 0 SEC 


. j — ] — L I — i I., j i i r— i i i i i i i i 

38 SB 70 98 1 18 130 150 170 198 210 230 258 278 290 318 330 350 370 390 ; 

MAGNITUDE OF ERROR. FT/SEC 



n m r 


MAGNITUDE OF ERROR. FT/SEC 


_L I- A i . I I . j I J I I I I I I 1 .1 

2H 30 36 42 40 54 60 66 73 79 05 91 97 103 109 115 121 

MAGNITUDE OF ERROR. m/SEC 


PILOT 17, a, = .15 SEC 


PILOT 18, o-, = .20 SEC 


I I I _L _L i _ I 1 — n — i -i i i i i i ii F I 1 I I I 1 I — I I i i i j i i i i— 

20 4 0 6 0 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 300 X 10 l ° 15 15 H5 75 105 135 165 195 225 255 205 315 345 375 405 435 465 435 J 

MAGNITUDE OF ERROR. FT/SEC MAGNITUDE OF ERROR. FT/SEC 

_l I I l J J I l J I I 1 J .1 J I I l I 1 I I I I I 

6 12 10 24 30 37 43 49 55 61 67 73 79 05 91 98 104 110 116 x 10' _ 5 5 1H 23 32 H 1 50 59 69 78 07 96 105 114 123 133 142 151 1 

MAGNITUDE OF ERROR. m /SEC MAGNITUDE' OF ERROR . m/SEC 


(e) Angle error definition a a = 0.20°. 
Figure B-4.- Concluded. 
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APPENDIX B - Continued 



MAGNITUDE OF ERROR, METERS 


MAGNITUDE OF ERROR, METERS 



MAGNITUDE OF ERROR. FEET 


I I . _L _L I 


-76 -70 - 64 - 58 -52 


~l J -I J . -L J l_ 1 . I _l l l l i i 

-40 -34 -27 -21 -15 -9 -3 3 9 15 21 27 34 40 46x10 s 

MAGNITUDE OF ERROR, METERS 


(a) Angle error definition cr a = 0°. 

Figure B-5-- Range error distributions at the end of the run (tg or ty) for 

the various pilots . 
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APPENDIX B — Continued 



MAGNITUDE OF ERROR, METERS 




_45^43 -40 -37 -34 -30~ - Z 7 ~^ 24 ~ l >\ -78 -15 -12 -9 -6 -3 0 3 6 9 12 I5xl0 5 

MAGNITUDE OF ERROR. METERS 


(b) Angle error definition a a = 0.05°. 
Figure B-5-- Continued. 
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(c) Angle error definition a a = 0.10°. 


Figure B-5-- Continued. 
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Continued 




MRGNiTunr or error, feti 
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1 -L _J 
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(d) Angle error definition cr a = 0.15°. 


Figure B-5-- Continued. 
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APPENDIX B — Continued 




MRGNITUDE OF ERROR. FEET 

I i i i i i i i i i 1 i i i 1 1 I ] I . j | 
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(e) Angle error definition 0 a = 0.20°. 
Figure Concluded. 
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